Competing interactions in two dimensional Coulomb systems: Surface charge 
heterogeneities in co-assembled cationic-anionic incompatible mixtures 
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A binary mixture of oppositely charged components confined to a plane such as cationic and 
anionic lipid bilayers may exhibit local segregation. The relative strength of the net short range 
interactions, which favors macroscopic segregation, and the long range electrostatic interactions, 
which favors mixing, determines the length scale of the finite size or microphase segregation. The 
free energy of the system can be examined analytically in two separate regimes, when considering 
small density fluctuations at high temperatures, and when considering the periodic ordering of the 
system at low temperatures (F. J. Solis and M. Olvera de la Cruz, J. Chem. Phys. 122, 054905 
(2000)). A simple Molecular Dynamics simulation of oppositely charged monomers, interacting 
with a short range Lennard Jones potential and confined to a two dimensional plane, is examined at 
different strengths of short and long range interactions. The system exhibits well-defined domains 
that can be characterized by their periodic length-scale as well as the orientational ordering of their 
interfaces. By adding salt, the ordering of the domains disappears and the mixture macroscopically 
phase segregates in agreement with analytical predictions. 



I. INTRODUCTION 

Biological and synthetic heterogeneous charged 
molecules are expected to self-organize in aqueous 
solutions into complex ionic structures. Co-assemblies 
of oppositely charged molecules are ubiquitous given 
that nucleic acids and most proteins are charged. 
The structure of oppositely charged biomolccular co- 
assemblies such as DNA-proteins in nucleosomes [|| and 
actin-protein complexes in the cytoskeleton 0], are the 
result of the competition of short range interactions, 
including excluded volume, and electrostatics. Cationic 
and anionic mixtures of lipids or peptide amphiphiles 
co-assembled into vesicles Q or cylindrical micelles 
[E IE Q are examples of co-assemblies stabilized by 
hydrophobic interactions and electrostatics. The sur- 
faces of such complexes of oppositely charged molecules 
may not be homogenous if the chemically co-assembled 
structures have net repulsive short range interactions 
among them, or if the charges exposed to to surfaces 
have different degrees of compatibilities with water. 
Understanding the surface assembly of a complex group 
of charged components may lead to a greater deal of 
understanding concerning the stability of self-assembled 
aggregates, or moreover, give insight into the complex 
behavior of lipid rafts and their contribution towards 
protein sorting and cell signaling |^. 

Bulk solution properties of electrostatic driven co- 
assemblies of cationic and anionic macroions have been 
extensively studied, such as DNA in cationic molecules of 
valence 3-1- and higher as well as other synthetic 

strongly charged polyelectrolytes in metallic multivalent 
salts llj. These hydrated multivalent ions are known to 
induce the precipitation of strongly charged chains of op- 
posite charge into dense ionic structures 12,13,14]. Co- 
assemblies of hydrophobic molecules of opposite charge, 
however, are less understood. Surface heterogeneities in 



co-assembled chemically incompatible oppositely charge 
molecules have been recently predicted analytically [Tjl- 
The surface charge heterogeneities are due to the com- 
petition between the long range electrostatic interactions 
(which decay as 1/r because the surface is embedded in 
a tree dimensional medium) and the short range inter- 
actions. The net incompatibility among the chemically 
different components of opposite charge promotes macro- 
scopic segregation. Electrostatic interactions, on the 
other hand, promotes mixing into an ionic crystal struc- 
ture. Consider cationic molecules with strong attrac- 
tions among themselves and restricted to surfaces such as 
cationic lipids adsorbed onto the surface of mica, which 
is negatively charged [Tsl [l6l| . The cationic molecules 
will aggregate into positively charged domains due to the 
strong net van der Waals attraction among them. The 
size of the domain, however, cannot grow past a charac- 
teristic size due to the high energetic penalty associated 
with the creation of a charged domain. This results in 
ordered finite size domains on the surface at low temper- 
atures 17]. Finite size charge heterogeneities have been 
observed experimentally on charged surfaces in the pres- 
ence of adsorbed self-aggregating molecules of oppositely 
charge Moreover, lattice Monte Carlo simulations 

of incompatible cationic and anionic molecules restricted 
to the surface of cylinders reveal many interesting finite 
temperature effects as well as various stripes structures 
along the cylinder at lower temperatures (iq . 

The formation of charged domains on a flat square lat- 
tice due to the competition between Coulomb interac- 
tions and net short range repulsion amongst oppositely 
charged molecules has been explored by simulation at 
zero temperature |19| ] and also by mean field arguments 
at high temperatures 0, 0] . These stoichiometric mix- 
tures develop ordered striped domains possessing a char- 
acteristic width that depends on the strength of the com- 
peting Coulomb and short range interactions at low tem- 
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peratures. At high temperatures percolated structures 
develop that resemble a spinodal decomposition pattern 
during phase segregation of binary systems, but growth 
is restricted, as in block copolymer systems with mi- 
crophase segregation [23, 0|. Here, we analyze the for- 
mation of the charged domains in two dimensions via 
molecular dynamics simulations at different ratios of the 
short and long range interactions. We analyze the sym- 
metric case of equal head size of stoichiometric mixtures 
of -1-1 and -1 charges with different effective interactions 
among them. Finite temperature effects are discussed. 
In Section II analytic arguments are given for the scaling 
of the surface charge domain size in various regimes of 
the degree of incompatibility. We justify the existence 
of charged domains in surfaces, as compared to a bulk 
three dimensional system. In section III we describe our 
simulations. In Section IV we discuss the results and the 
last section we give the conclusions. 



II. THEORY 

The phase behavior of the ionic mixture can be exam- 
ined in a simplistic manner analytically in two separate 
regimes. At higher temperatures, we consider small den- 
sity fluctuations around the mean. At low temperatures, 
when the system exhibits strongly segregated domains, 
we assume the system is periodic. At low temperature 
values, or high values of the magnitude of short range at- 
traction, the system exhibits well-defined periodic lamel- 
lar when the charge surface coverage of the positive and 
negative molecules are equal and confined to a flat sur- 
face. The free energy of the system is dominated by the 
electrostatic cohesive energy in addition to the interfacial 
contribution to the free energy, which is characterized by 
the line tension, 7 per thermal energy fc^T. Within the 
strong segregation regime, the entropic contribution to 
the free energy is negligible. 

Following the example of the free energy for a incom- 
pressible two dimensional system of a mixture of positive 
and negative components |l7l |. we generalize the results 
for a course grained free energy scaling analysis for a d 
dimensional system of Na positively and Nb negatively 
charged components interacting with a three dimensional 
Coulombic 1/r potential. The free energy can be writ- 
ten as sum of the total electrostatic interactions and the 
contribution from the line tension of each periodic segre- 
gated domain. Each charged domain is approximated by 
a electroneutral unit cell which has a characteristic lat- 
tice length L, dimensions L"^, and an associated charge 
density a. The net free energy per total number of par- 
ticles N = Na + Nb, in units of kBT, can be written 
as 

'cell a'' f lBa\s2{L'')^\ fFoa'' 
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Here, si and S2 are geometrical parameters that depend 
on the characteristic geometry of the underlying unit cell. 



a'^ represents the size of the particle, and Nceii represents 
the number of particles per unit cell. The Bjerrum length 
Ib is given by. 
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Fo and Lo are system dependent parameters, defined by 
the minimization of the free energy of the system with 
respect to the characteristic size of the system, L, 



Fa^ 
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and 
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The dimensionless free energy per unit area in terms of 
si and S2 is then described by 



(5) 



where D ~ L / is the ratio of the characteristic size of 
the unit lattice to the length of the system. Minimizing 
the dimensionless free energy with respect to D gives the 
free energy of the favored periodic structure as 
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Depending on the area fraction of charge coverage, /, 
and the geometry of the unit lattice cell, the free energy 
can be calculated for different sets of crystalline struc- 
tures. For an ideally symmetric system, consisting of 
equal components of positively and negatively charged 
molecules with similar head group sizes, / is 1/2. The 
minimum free energy in this case, for a two dimensional 
system, is characterized by lamellar structures. 

We consider a line tension that is proportional to the 
immiscibility of the component molecules, x ■ The Flory- 
Huggins paramater, x, is defined as the difference in the 
magnitudes of the short range interactions between two 
components as x = (ei2 — 5(^11 + ^22)) /^bT, where the 
Cij represents the pair interaction energy between i and 
J. For a lower, or two dimensional system, Lq, would 
be comparably larger than for a three dimensional sys- 
tem due not only to the l/d power law dependence but 
also to the decreased value of the Bjerrum length Ib for 
BpSurface in contact with water. For a surface in con- 
tact with an aqueous solution the mean permittivity of 
the medium is much higher than in a dense three dimen- 
sional system, which decreases the Bjerrum length Ib, 
and thus the magnitude of Lo, significantly. For these 



3 



reasons, patterning on a surface due to the competition 
of electrostatic interactions with short range interactions, 
is considerably more feasible than the creation of charge 
domains in a bulk three dimensional system. 

Comparing length scales with experimental systems, 
consider a two dimensional system of a single layer of 
positively and negatively charged lipids at an interface 
between water and an alternate medium. The average 
dielectric permittivity of at the interface ei ^ AO, in be- 
tween that of the water t^ater ~ 80 and that of the dense 
medium tmedium ~ 1- This would correspond to a Bjer- 
rum length Ib ^ 2nm in terms of the a classical elec- 
trostatic interaction between charged head groups of the 
lipids exposed to the aqueous interface. Considering a 
large magnitude of the net interaction between tails of 
interacting lipids at the interface (x ~ 15), depending 
on the length of the hydrophobic tail of the molecules 
20 carbons) and the charge density of the head-group 
.6/nm^), this could correspond to a fairly large equi- 
librium domain size Lo (~ 80nm). Domains of this size 
or larger have been seen for experimental systems of com- 
peting short range and long range electrostatic interac- 
tions, although in comparing with these systems, a vari- 
ety of kinetic and specific interaction effects should also 
be considered [2ll2^. 

Next, consider the opposite, high temperature regime. 
Since the system does not exhibit well-defined periodic 
structures, the entropic contribution to the free energy 
cannot be ignored. In this case, linear response theory 
or the Random Phase Approximation for a compressible 
binary systems 0| is used to describe the behavior of 
the correlations as a function of the relative strength of 
the short range attraction and the electrostatic interac- 
tions. The Random Phase Approximation involves an 
expansion of the free energy of the system in terms of 
density fluctuations, neglecting all terms of larger than 
second order. For a general system of N components, 
where i and j represent components of a different type, 
the partition function can be written as jlll l25j | 



Hiickel potential. 
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where the Hamiltonian of the system is represented by 



(1) 
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where k, the inverse screening length, is defined by the 
concentration of salt in the solution. We assume that the 
density is a smooth function and can be represented by 
the sum of its Fourier components 



In this case, the partition function becomes 



(12) 
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where A represents the area of a two dimensional plane 
in a three dimensional volume V. Zq includes the k zero 
and the self energy terms. U^"' is the sum of the in- 
teraction energies of the system, consisting of the short 
range interactions due to the excluded volume and hy- 
drophobic interactions, vfj^{k), as well as the long range 
electrostatic potential, u|j(fc). 

For an incompressible system of i same-sized compo- 
nents, 



(14) 



For the case of a incompressible, neutral, symmetric sys- 
tem we also assume that p+{k) = —p-{k). The electro- 
static potential is the two dimensional Fourier transform 
of the screened Coulomb interaction between charge den- 
sity fluctuations. 



Uel{k) = 
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where a represents the charge density of the system and 
zt represents the total positive and negative charge of the 
components. In this case, the inverse structure factor has 
the following contributions. 



It is assumed that the interparticle potential can be bro- 
ken up into a short range and long range electrostatic 



potential, Vi. 



,SR 



The short range contribution 



is assumed to be of the form of the Fourier transform of a 
Gaussian potential, which has been shown to reason ably 
predict thermodynamic properties of binary systems [26j , 



<"M = ^e--/-. (10) 
The long range potential is represented by the Debye 
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^o(fc) 
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-2x + xk^ + Uei{k%16) 



The structure function has a peak at the most probable 
wave lengths, k* . For the case when there is no screen- 
ing the location of the peak, k*, scales with the Bjerrum 
length, Ib and magnitude of short range attraction, e, 
as k* ^ (e/Zs)'^. The scaling of the periodic order of 
the system changes at the transition temperature from 
— (1 + d) at higher temperatures considering small den- 
sity fluctuations to —d at lower temperatures (Eq. 4), 
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which is predicted using the previously described the- 
ory of strong segregation. For a two dimensional system, 
which is the subject of interest, the scaling is predicted to 
change from —1/3 to —1/2 as the temperature decreases. 

At high temperatures, in the nearly isotropic state, the 
total free energy of the system per unit volume, in units 
of fcsT, can be written as 



TABLE I: Interparticle Potentials 



Interactions 




+ 
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Uc- 


-ELC + UlJ 


Uc-ELC + UhC 




Uc- 


-ELC + UhC 


Uc-ELC + UlJ 



AkgT 



01n(/) + (1 - 0) In (1 - 0) - + F,i,/{kBTlX7) 



where A represents the two dimensional area of a plane 
and Ff,ie/ {ksT) represents the one loop corrections ob- 
tained by integrating the charge density fluctuations [2^ . 



III. MODEL AND SIMULATION DETAILS 

The model system is composed of a mix of N+ pos- 
itively and iV_ negatively charged monomer units in a 
simulation box of size Li^ . The molecules are confined 
to a two-dimensional plane perpendicular to the Z axis, 
with periodic boundary conditions in the X and Y direc- 
tions. Each monomer unit represents a charged biological 
or polymeric molecule, that interacts attractively with a 
like monomer via hydrophobic forces. In this paper, only 
symmetric mixtures are considered, where the charge and 
radius of the positively and negatively charged monomer 
units are equivalent. The total system is electroneutral. 
We are interested in the case where the two-dimensional 
layer exhibits well-defined periodic patterns along the 
surface of the plane. Fluctuations perpendicular to the 
monomer plane are restricted. 

Constant N,V,T Molecular Dynamics simulations were 
performed using Espresso, simulation code developed by 
the MPIP-Mainz group of Polymer Theory and Simu- 
lation (http://www.espresso.mpg.de). A stochastic or 
Langevin thermostat is used, to ensure a constant tem- 
perature, along with a Verlet algorithm to calculate par- 
ticle velocities at each timestep. The unit of energy is 
e, of length cr, and of mass m. Temperature is then de- 
fined in terms of e/fc^T and time in units of y^r^m/e. A 
full Coulomb potential is used for calculations of charge- 
charge interactions. The ELC (Electostatic Layer Cor- 
rection) method developed by Arnold et al. to sum 
the electrostatic energy contribution to the free energy 
[27i . This method is a correction to the P3M Ewald 
summation technique "29] , in which the Fourier transform 
of the electrostatic contribution to the energy is summed 
using a mesh formulation. In addition to full electrostat- 
ics, a case of a screened Debye-Hiickcl interaction is con- 
sidered to look at the subsequent melting of the periodic 
structures when the potential is screened. Table 1 sum- 
marizes the interaction potentials between the positive 
and negative component monomers in the system. The 
potential between charges is the full Coulomb potential, 



Uc- 



ELC 



lBTqiq2 



where 1^ represents the Bjerrum length of the system. 
For the present simulation results, Ib's of O.lcr, 0.2a and 
0.5ct are considered. Considering an average dielectric 
permittivity of the medium (e^ ~ 80) this corresponds to 
a fairly large headgroup size {^-^ 20 A). The short range 
interaction between like monomers is the classic Lennard 
Jones potential, 



Ulj = 4e 



(7) 
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r <rr. 



(19) 



where a is the monomer radius, and the potential is cut 
at a radius Tc of 2.5a. An additional term is also added 
to the potential to keep the derivative continuous at Tc. 
Uhc is the same as Ulj, with a cutoff radius, of 2^/^a, 
including only the repulsive part of the potential, which 
represents the excluded volume of the molecule. 

Initially, a fairly dense surface density, p, of 0.6 was 
considered, to compare with phase behavior predicted 
by strong segregation theory, while remaining sufficiently 
far from the two dimensional hard disc crystallization 
regime of approximately p = 0.89 determined by previous 
MC and MD simulations 's^. This also allows sufficient 
diffusion for the system to equilibrate. 



P 



{N++N^)TTa^ 
4L2 



(20) 
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Phase behavior in comparison with theory at lower sur- 
face densities is slightly more complex and will be dis- 
cussed in a later paper. The majority of simulation re- 
sults are presented for a system size of 1000 charged 
monomers, while finite size effects are explored by in- 
creasing the system size by a factor of 2. Approximately 
10^ MD steps are used to equilibrate the system; the equi- 
libration time grows increasingly longer at higher values 
of the Bjerrum length. 



IV. DISCUSSION OF TWO DIMENSIONAL 
PHASE BEHAVIOR 

At lower values of e, domains of positive and negative 
component monomers appear in the system. As the mag- 
nitude of e increases, the domains begin to increase in size 
in an isotropic manner, forming a percolated structure. 
As the value of e further increases, the domains begin 
elongate and then to orient into well-defined lamellar, 
breaking the symmetry of the system. Increasing even 
further, the lamellar begin to widen. Average internal 
energy and heat capacity per particle are calculated at 
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FIG. 1: Snapshots of the system at e = 

1.0(a), 2.5(6), 4.0(c)fcBT at a constant Ib of 0.2a. In- 
troduction of K = 20a (d) induces macroscopic phase 
segregation at an e of 4.0fcsT, Ib of 0.2a. 



several values of Bjerrum length {Ib — O.Ict, 0.2(t, 0.5a). 
At higher values oi Ib, electrostatics plays a more im- 
portant role in the equilibrium configuration of the sys- 
tem. The electrostatic repulsion between like charged 
monomers increases. In order to minimize this contribu- 
tion to the free energy, the stripes become thinner. The 
average internal energy {< E > /N) and heat capacity 
per particle (Cv) are calculated at two different values 
of the Bjerrum length {Ib = 0.1a, 0.2cr, 0.5cr). The aver- 
age internal energy is less negative at the higher value of 
Ib due to the increased repulsion between like charged 
head groups. At lower values of Ib the heat capacity dis- 
plays a peak, which corresponds to the crossover from the 
percolated phase to the lamellar phase. The magnitude 
of this peak increases and shifts to the left as the value 
of the Bjerrum length is decreased. Lamellar spacing is 
systematically characterized by the calculation of the two 
dimensional structure factor, S{k), where r corresponds 
to a vector in the x,y plane. 

5(fc)= y"5(r-r^)e*'=-V^-'''dV (21) 

S{k) displays a peak at k* corresponding to the inverse 
lamellar spacing in the direction perpendicular to the 
lamellar. As a function of e, the peak location corre- 
sponds to scaling predictions by strong segregation the- 
ory at high values of e (fc* ^ e~^/^). At lower values, the 
location is consistent with predictions by the Random 
Phase Approximation (fc* ~ e^^^'^). The orientational 
order of the domains can be characterized by the inter- 
facial orientational order parameter g2 [sij , 

i=i j=i 



where Ni is the number of neighbors of opposite type 
of monomer at an interface and Oij is the angle between 
two neighbors. A neighbor is defined as two particles 
of different type, within range of short range attraction 
{rij < Tc). As the magnitude of short range attraction in- 
creases, the calculated order parameter increases, which 
corresponds to the ordering of the domains by the devel- 
opment of orientational order at the interface. The in- 
crease in order of the domains, indicated by an increase in 
the order parameter 52 1 proceeds the location of the peak 
in the heat capacity. Higher values of the Bjerrum length 
Ib correspond to a higher value of the order parameter 32 
for stronger short range attraction. As the electrostatic 
contribution to the segregation increases, the character- 
istic domain size decreases, but the orientational order of 
the domains increases. At lower values of Ib, this initial 
increase is followed by a levelling off, or slight decrease, 
that corresponds to the formation of holes at the inter- 
face. The holes disrupt the orientational order, or the 
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FIG. 2: (a) Heat capacity per particle (Cv) and (b) average 
internal energy {< E > /N) at several values of the Bjerrum 
length {Ib = O.lcr, 0.2(7, O.Scr). The heat capacity displays 
a peak, which corresponds to a crossover from percolated, 
random domains to a lamellar phase. The magnitude of the 
peak increases and shifts to the left as the value of the Ib is 
decreased. 
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FIG. 3: (a) The interfacial orientational order parameter gi 
at several values of the Bjerrum length {Ib = O.lcr, 0.2(t, 0.5cr) 
as a function of e. As the orientational order of the domains 
increases, (72 increases from to a finite value, (b) The location 
of the peak k* in the structure factor S(k) as a function of e. 
The scaling of fc* with e changes from —1/3 to —1/2. 



is considered by using a screened Debye Hiickel potential 
instead of the Coulomb potential for electrostatic inter- 
actions, 



DH 



lBTqiq2e' 



(23) 



where k represents the screening length due to the sur- 
rounding three dimensional solution of ions (see Equa- 
tion 10). At higher values of electrostatic screening 
(k = 5cr, 10(7, 15(t), examining the behavior of the sys- 
tem with an intermediate value of short range attrac- 
tion (e = 2.0), the system phase segregates into two 
macroscopic charged domains of positive and negative 
ions. The peak in the structure factor indicates that the 
segregation length-scale is nearly constant as a function 
of the screening length. This is in agreement with ana- 
lytic theory JT?] . The location of the peak shifts to lower 
values with an increase in the characteristic size of the 
simulation box. To determine if these simulation results 
are consistent with theoretical predictions, we examine 
the behavior of the inverse structure factor (Eq. 16) in 
a regime where the scaling of the peak in the structure 
factor from simulation results is still consistent with lin- 
ear response theory. We find that as k, the magnitude of 
screening by the ions of solution, increases, the value of 
q* goes continuously to zero. 



X 



1/2 



(24) 



before the structure factor diverges, at which there is 
macroscopic phase segregation. At higher values of short 
range attraction, the structure factor diverges when q* > 
0. This is in agreement with analytical predictions from 



hexagonal packing, of the monomers within the segre- 
gated domains. This is equivalent to the inclusion of a 
ternary component with a non-selective interaction be- 
tween positive and negative components. Initial exam- 
inations on the finite size effects of the system are ex- 
plored to determine the effect of the periodic boundary 
conditions on the ordering of the more strongly segre- 
gated lamellar. Doubling the size of the system at larger 
values of the short range attraction (e = 4A:bT), quanti- 
tatively affects the ordering of the lamellar by decreasing 
the alignment of the domains along the boundaries of 
the system and increasing the fluctuations along the in- 
terface. This results in a characteristic decrease in the 
order parameter, 32 from .43(±.02) to .37(±.02). Fur- 
ther system sizes were not examined due to sufficiently 
high surface density; the calculation of the electrostatic 
energy is slow to converge at these density ranges. 

Introduction of electrostatic screening, or including the 
effects of high salt on the local ordering of the system. 




FIG. 4: The location of the peak, k* , in the structure 
factor, S{k), as a function of k as predicted by linear re- 
sponse theory for several values of short range attraction 
(x = 0.1,0.5,2.0,5.0,20.0) for an intermediate strength of 
the electrostatics, Ib = 0.2. The shaded area indicates S{k) 
diverges at a finite value of k. 
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FIG. 5: The location of the peak, k* , in the structure factor, 
S{k), as a function of k from simulation results spanning a 
range of the screening parameter {k — 3 ~ 15a) at an inter- 
mediate value of short range attraction (e = 2.0). Increasing 
the size of the system from L = 36 to L = 51 decreases the 
average value of k* . 



the strong segregation regime, which predicts a discon- 
tinuous jump from finite sized periodic cells to an infinite 
cell at a value of k which is inver sely proportional to pe- 
riodic length-scale of the system |l7l |. 



V. CONCLUSIONS 

Molecular Dynamics simulations of oppositely charged 
monomers, interacting with a short range LJ potential 
and confined to a two dimensional plane, are examined 
at different strengths of short range attraction and long 
range electrostatics. The system exhibits well-defined do- 
mains; the width and ordering of the domains are depen- 
dent on the depth of the LJ well, e, and the strength of 
the Coulomb interactions, Ib- The length-scale of the or- 
dering of the system can be quantitatively characterized 
by the two dimensional Fourier transform of the density, 
S(k), where k is the inverse spacing of the system. S{k) 
has a peak k* which scales with the line tension of the 
domains, 7. The underlying assumption of strong seg- 



regation theory is that the microphase regions of charge 
are well defined and periodic, with a line tension 7 that is 
proportional to x- Since the magnitude of the short range 
attraction e is proportional the Flory-Huggins interaction 
parameter x, k* should scale with e in the regimes where 
strong segregation theory holds JL7|. It is shown that, 
at higher values of e, the scaling of k* with e is consis- 
tent with theory. At lower values of e, a different scaling 
is found, which is consistent with that which is found 
using linear response theory. Electrostatics represents a 
more important contribution to the characterization of 
the interfacial line tension in this regime. 

The degree of ordering can be examined by the calcu- 
lation of the interfacial orientational order parameter, 172 ■ 
The transition from a random, percolated domain struc- 
ture to well defined lamellar is a gradual transition, that 
is demonstrated by the gradual increase of the parameter 
(72 as a function of e. This result is consistent with what 
one would exhibit with a Kosterlitz and Thouless type 
transition 30], in which the two dimensional ordering 
the system exhibits a continuous phase transition, that 
can be defined by a similar positional order parameter. 
Initial examinations of the finite size effects of the sys- 
tem indicate that the degree of ordering is slightly influ- 
enced by the periodicity of the simulation box, however, 
further examinations were not made due to the compu- 
tational intensiveness of the electrostatic energy term. 
Decreasing the strength of the electrostatics in the sys- 
tem, by changing the charged interaction from a straight 
Coulomb potential to a screened Debye Hiickel interac- 
tion, the ordering of the system disappears and the mix- 
ture phase segregates, which is consistent with analytical 
arguments. 
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